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ABSTRACT 


A quadrilateral vortex-lattice method was applied to a formation of three 
wings to calculate force and moment data for use in estimating potential bene- 
fits of flying aircraft in formation on extended range missions , and of antici- 
pating the control problems which may exist* The investigation led to two types 
of formation having virtually the same overall benefits for the formation as a 
whole, i.e. f a 11 V” or echelon formation and a double row formation (with two 
staggered rows of aircraft) * These formations have unequal savings on air- 
craft within the formation, but this allows large longitudinal spacings be- 
tween aircraft which is preferable to the small spacing required in formations 
having equal benefits for all aircraft* A reasonable trade-off between a prac- 
tical formation size and range benefit seems to lie at about three to five air- 
craft with corresponding maximum potential range increases of about 46% to 67% * 
At this time it is not known what fraction of this potential range increase is 
achievable in practice. 
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FORMATION FLYING BENEFITS 

BASED ON VORTEX-LATTICE CALCULATIONS 

* 

Brian Maskew 


SUMMARY 


A quadrilateral vortex- lattice method was applied to a formation of three 
wings to calculate the force and moment data for use m estimating potential 
benefits of flying aircraft m formation on extended range missions, and of 
anticipating the control problems which may exists 

The investigation led to two types of formation having virtually the same 
overall benefits for the formation as a whole. Calculations indicate that "V" 
or echelon formations have an induced drag reduction of 80% on all but the 
leading aircraft. Double- row formations (with two staggered rows of aircraft) 
have virtually no induced drag savings m the leading row, but aircraft in the 
second row have an induced thrust of 48% of their free-air induced drag. Both 
formations have large longitudinal distances between aircraft (about three wing 
spans), and are, therefore, preferable to formations having equal benefits for 
all aircraft. The latter formations require vary small spacings as observed in 
migrating bird formations. The double-row formation has an advantage over the 
M V ,f formation in that rolling moment required for trim is zero. Both forma- 
tions are insensitive to longitudinal spacing, but benefits (and trimming prob- 
lems) decrease rapidly with lateral and vertical movements from the optimum 
locations . 

A reasonable trade-off between a practical formation size and range bene- 
fit seems to lie at about three to five aircraft with corresponding maximum 
potential range increases of about 46% to 67%. At this time it is not known 
what fraction of this potential range increase is achievable m practice. 

The optimum locations require the trailing vortices from the leading wings 
to pass close to the tips of the followers, and this poses a problem for exist- 
ing calculation techniques. The problem affects not only the prediction of the 
maximum benefits, but also the stability and control aspects of holding forma- 
tion. Further work, therefore, is recommended to investigate force and moment 
calculations with close vortex/wing interactions and to refine the calculation 
method. The results should then be used to analyse the control, stability and 
range performance of specific aircraft in formation. 


* 

Senior Research Scientist, Analytical Methods, Inc., Bellevue, Washington 
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INTRODUCTION 


Recently there has been renewed interest in the concept of flying aircraft 
in formation to increase their range or to reduce the amount of fuel required 
for a given range. Modern avionics and control augmentation systems (including 
alternative autopilot modes) would facilitate holding accurate formation posi- 
tions over long periods and would thereby help to achieve the maximum benefits. 
The benefits arise from an induced drag saving that is caused by favorable 
aerodynamic interference between the aircraft* So far, formation flying calcu- 
lations, such as references 1 through 4, have been based mainly on a simple 
horseshoe vortex representation of each wing; Shollenberger's calculations 
(ref. 4) , however, indicate important differences between the results from that 
simple approach and from one based on elliptic loading when the spanwise spac- 
ing between wing tips is less than about a semispan. In fact, significant 
benefits are obtained only with spanwise spacmgs considerably smaller than a 
wing semispan, and so a detailed wing representation is clearly required for 
accurate calculations. This requirement is especially important in practical 
calculations in which each aircraft is trimmed in roll. Formations considered 
in previous calculations were mainly based on achieving the same induced drag 
saving for each aircraft m the formation — a situation that evidently exists 
in migrating bird formations. For aircraft, however, such formations require 
dangerously small spacmgs. 

The present objective is to apply a vortex- lattice method (ref. 5) to cal- 
culate formation flying benefits and to consider a range of formations to find 
safer alternatives to the "equal benefits" formation. Use of the vortex-lat- 
tice method allows a more realistic calculation than is possible with the 
simple horseshoe- vortex model? the mutual interference effects can be represen- 
ted in greater detail, and each wing can be trimmed in roll. THe~mam aims of 
this study are to find optimum formation configurations to assess the loss in 
benefits caused by small movements from the optima, and to assess the degree 
of sensing and control required to maintain range benefits within given limits. 
The wmg planform used in the calculation is based on a possible military 
S.T.O.L. transport aircraft configuration. 
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LIST OF SYMBOLS 


x, y, z Cartesian coordinate system (see fig. 1) 

C lift coefficient ~ lift/(g s) 

Li CD 

C D induced drag coefficient = induced drag/^q^s) 

i 

rolling moment coefficient = rolling moment/ (q^sb) 
free- stream dynamic pressure 


S 

a 

€ 

b 


true air speed of formation, m/s 
wing planform area, m 2 

wing incidence relative to the flight vector, deg. 
flaperon deflection, positive downwards, deg. 
wing span, m 


Subscripts: 

A inner wing 

B outer wing 

S single 

F formation 


See figure 1 
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CALCULATION METHOD 


Basic Method 


The method (ref. 5) used here is based on vortex-lattice theory (e.g., 
ref. 6), but the wing lattice is formed into quadrilateral vortices (fig. 1) 
instead of horseshoe vortices. Each quadrilateral has a control point at which 
the boundary condition of zero normal velocity is specified; the normal veloci- 
ty includes a component of the free- stream velocity and contributions from all 
the trailing vortices and quadrilateral vortices in the lattices (i.e., on all 
the wings m the formation in this case) . This gives a set of simultaneous 
equations in the unknown quadrilateral-vortex strengths; for example. 


1 _ 
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where 


,, , th 

is the k 


quadrilateral- vortex strength. 


A. is the influence coefficient for the normal component of 
^ velocity (ref. 5) induced at the j*** 1 quadrilateral control 
point by the kth quadrilateral vortex; 


B 

3 


is the normal component of free -stream velocity at the 3 
control point; and 


N is the total number of quadrilateral vortices in the system* 


Initially, the trailing vortices are assumed semi-infinite in the stream- 
wise direction, but the method incorporates an iterative procedure for calcu- 
lating the trajectories of the vortices for a force- free wake* In this pro- 
cedure r the first part of each trailing vortex is divided into a number of 
short, straight segments, and each segment is aligned with the local mean ve- 
locity vector* As the wake geometry changes, the trailing vortex contributions 
to the coefficients. A, in Equation (1) change, and so a new vortex strength 
solution is calculated at each iteration. 

For most of the present calculations, the procedure was terminated at the 
end of the first pass, i.e*, with the straight undeflected wake. One case was 
attempted with the wake iteration to assess the effect of wake roll-up. 

When the quadrilateral vortex strengths are known, the forces and moments 
are calculated by applying the Kutta Joukowski law to each bound vortex segment 
m the lattices, viz.: 


2V sf 


( 2 ) 



PLANE OF SYMMETRY 



Figure 1. Vortex-lattice Arrangement and 


NOTE: z IS POSITIVE 

UPWARDS 


NTROL POINT 


c=- y 



Definition 




where 


Cg is the force vector coefficient which is resolved into the 
lift and drag direction; 

V is the calculated local velocity vector at the bound vortex 
segment mid-point; and 

£ is the vortex segment vector. 


The force is assumed to act at the segment mid-point for the purpose of 
evaluating moments. 


Formation Representation 


Symmetrical formations are assumed; this halves the total number of un- 
known quadrilateral vortex strengths for a given problem. The formations con- 
sist of three wings and each wing has a lattice of 4 (chordwise) by 18 (span- 
wise) quadrilaterals, figure 1. The total number of unknowns, therefore, is 
108. 


Ideally, more quadrilateral elements should be used, e.g., 6 x 30 on each 
wing, but the computer program is limited to 150 unknowns and 30 trailing vor- 
tices at this time. Also, the lattice should be set in from each wing tip by 
a quarter of a quadrilateral "span" for best results (ref. 5) ; in the program 
coding prevented this from being applied directly on the inner tips of the 
outer wings, and so tip inset was not used in these calculations. Because the 
main interest in this study is m ratios of induced drag, etc., (i.e., forma- 
tion values/free-air values) rather than absolute values, neither of the above 
factors will be a serious drawback to the present calculations. 


Trimming Procedure 


Each wing is trimmed m pitch so that all the wings in the formation have 
the required lift coefficient, C^. The outer wings are also trimmed in roll 
using a downward-deflecting "flaperon” (fig. 1) . It is assumed that the flap 
mechanism could be arranged to allow the flaps on one side only to be deflected 
to provide the basic asymmetry required in, some of the formations; the aileron 
would be treated similarly, and would be used for periodic trimming purposes 
to maintain formation position. 

In the calculations, the wings are set initially at incidence a (i.e., 
a A = otg = a) giving the required lift coefficient in free-air conditions, and 
the outer wings have zero flaperon deflection, i.e., = 0. The initial solu- 
tion gives lift coefficients and Cjjg for the inner and outer wings, res- 

pectively, and an out-of -balance rolling moment, C^ B , for the outer wings. 

Using superscripts o and n for the old and new conditions, respectively, the 
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The trimming steps are repeated until the value and zero rolling moment are 
achieved. u 
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BASIC DATA 


Wing Geometry 


Figure 1 shows the wing planform including the flaperon. The basic de- 
tails are given in table 1. For one calculation the flaperon span was extended 
over the full semispan* 


Formations 


The coordinate system for defining the formations is shown in figure 1, 
and figure 2 shows a summary of the formations considered* Wing A, the inner 
wing, is fixed at the origin, while Wing B (and its image) moves along five 
scan lines, one streamwise, two spanwise and two vertical* The streamwise scan 
is in the plane z = 0, and goes from three spans upstream of the origin to 
three spans downstream. Except for one case where all three wings are tip to 
tip, the spanwise position of Wing B in the streamwise scan is at y/b = 0.89; 
this places the tip of the following wing near to the centroid of shed vor- 
ticity from the leading wing semispan. 

There are two spanwise scans m the plane z - 0, one at three spans up- 
stream of the origin and one at three spans downstream. Spanwise positions 
range from y/b = 0.83 to 2 m each scan. 

There are two vertical scans, one at three spans upstream and the other 
at three spans downstream of the origin. These scans include the common 
points in the streamwise and spanwise scans, viz., y/b = 0.89. The vertical 
movement from z = 0 is 0.5b. 

In formations where the trailing vortex- lattice from the leading wing 
overlaps the lattice on the following wing, (i.e., formations with y/b < 1) 
a numerical problem occurs in the calculations if the lattices do not match 
(ref. 7) ; only spanwise positions which line up the lattices are therefore 
considered (see fig. 1) . 


Single Wing Aerodynamic Characteristics 


Figure 3 shows the basic lift and induced drag characteristics calculated 
for a single wing represented by the 4 x 18 lattice. Because of the limita- 
tions mentioned under "Formation Representation", these characteristics should 
be regarded as approximate m the absolute sense; they are used for trimming 
the calculations here and for the base m the ratio between formation and free- 
air conditions. (However, a calculation was made at 0t = 7° with a tip inset of 
a quarter quadrilateral span and a 4 x 40 lattice, i.e., a 4 x 20 lattice on 
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SCAN LINES 
(LOCUS OF POSITIONS 
FOR WING B) 
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i 

\ 



POSITIONS OF 
VERTICAL SCAN LINES 




Figure 2. 


Summary of Formation Positions Considered 



Table 1. Basic Wing Geometry 

(Based on a Semispan of 1.0) 


span 2.0 

geometric mean chord 0.25 

area 0.5 

aspect ratio 8.0 

taper ratio 0.33 

sweepback (quarter- chord line) ... 5° 

dihedral 0 

twist 0 

flaperon span . . 0.44 to 1.0 

flaperon chord/local chord 0.19 
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the half span. The results (fig. 3) show only a small displacement from the 
4 x 18 lattice characteristics.) 

The lift curve slope to be used in equations (3) and (5) is, therefore: 


dC L 

-r— = 0.0851 per degree (6) 


The f laperon effectiveness (for equations (4) and (5) ) calculated at a = 0 
with the 4 x 18 lattice is: 


= 0.002982 per degree 

(7) 


= 0.01060 per degree 


(The corresponding quantities calculated with a 4 x 20 lattice were 0.003014 
and 0.01065, respectively. Because of asymmetry, this calculation considered 
the full wing and the lattice was, therefore, limited by the number of trail- 
ing vortices allowed m the computer program, viz., 30 at this time.) 


dC & 

d§ 


d 5 


Data Points 


Calculations performed for each formation configuration (after trimming) 

evaluate the induced drag for each wing and for the formation, together with 

the incidence (for the required lift coefficient) and the f laperon angle. The 

induced drag and the incidence are then divided by the corresponding free-air 

values (at the same C_ ) . The single-wing data are : 

Xj 

for C = 0.6, 

X; 

C D = 0.01296 

l 

a = 7.051° 
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and for C = 1.2, 

L 

C D = 0.05189 
i 


a = 13.92° 


(Both sets of data are for the 4 x 18 lattice.) 

A few data cases considered at the higher C L values gave essentially the 
same incidence and induced drag ratios as for the C L = 0.6 case because of the 
approximate linear relationships of C ~ a and C ~ C 2 (fig. 3) . 

Xj D . L • 

The incidence and flaperon angles for each formation are related back to 
approximate lift and rolling moment increments, respectively (using the deriva- 
tives from equations (6) and (7)). These quantities vary approximately linear- 
ly with Cjj and are presented divided by the C L value of the calculation (i.e., 
0 ., 6 ) . 
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RESULTS AND DISCUSSION 


Streamwise Scan 


Figure 4(a) shows the induced drag ratio variations as wing B (and its 
image) moves from three spans upstream to three spans downstream of wing A. 

The corresponding variations in incidence ratio, lift increment, flaperon angle 
and rolling moment are given m parts (b) , (c) (d) and (e) of figure 4, res- 
pectively. 

The position for equal induced drag benefit for all three wings occurs at 
a very small x/b value, i.e., 0.15, giving a streamwise distance of approxi- 
ma tely 1_5 tip chords between the tip trailing-edge of wing A and the tip 
leading-edge of wing B. 

The induced drag ratio for the formation of three wings remains virtually 
constant at 0.468 over the complete scan (and therefore complies with flunk' s 
stagger theorem) . The individual wing induced drag values change markedly in 
the region where wings A and B change place as leader, but, by three spans 
upstream and downstream of the origin, the values have almost reached steady 
conditions. Wing B achieves an induced drag ratio of 0.21 at x/b = 3.0; the 
limiting value is 0.20 for the larger downstream distances assuming the total 
for the three wings remains constant and that wing A goes to free-air con- 
ditions . 

At x/b = -3, wing B has an induced drag ratio of 0.94, but wing A is now 
m an upwash field from both wings B and receives a negative induced drag ratio 
(i.e., a thrust) of -0.48. 

The angle of incidence for a given C L is reduced appreciably in formation, 
figure 4 (b) . The incidence ratio for wing B is 0.83 at x/b = 3 while that for 
wing A with x/b - -3 goes down to 0.767. The corresponding lift increment 
ratios (before trimming) are 0.127 and 0.235 respectively (fig. 4(c)). 

The flaperon angle required to trim the rolling moment on wing B is 2.7° 
when x/b = 3 (fig. 4(d)), and represents an induced rolling moment coefficient 
of -0.01330^ (fig. 4(e)). Formations with wings B ahead of wing A have almost 
negligible rolling moment. 

One case was considered at x/b = 3 with a full semispan flaperon to see 
if this would reduce the roll trim drag. The larger flaperon, rather than 
reducing trim drag, gives a small increase (fig. 4(a)). (But the smaller de- 
flection required by the larger flaperon should reduce any profile drag 
penalty.) Evidently, deflecting the shorter flaperon must have improved tne 
basic loading distribution towards the ideal elliptical shape, giving a 
negative trim drag for this wing. Deflecting the larger flaperon essentially 
restored the basic loading distribution and induced drag. 
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(a) Induced Drag Ratio. 


Figure 4. 


Effect of Streamwise Position: z/b = 0.0; 

y/b = 0.89 (except y/b = 1.0 at x/b = 0.0). 



Figure 4. Continued. 








Spanwise Scans 


The induced drag ratio variation with spanwise movement at x/b = 3 is 
shown in figure 5(a) for wings A and B and for the formation of three wings* 

The minimum value of 0.21 for wing B at the spanwise position of 0*89 (fig. 
4(a)) increases to 0.28 at y/b = 1.0, then increases rapidly and is 0.49 at 
y/b = 1*05* The variations in the incidence, induced lift, flaperon angle and 
rolling moment ratios, figure 5(b), (c) , (d) and (e) , respectively, become very 

steep as y/b decreases towards 1? i.e., as the inner tip of wing B approaches 
the edge of wing A's trailing lattice. As y/b increases from 1.0 to 1.05, the 
incidence ratio increases by 0.07 (fig. 5(b)), the induced lift decreases by 
0.0370^ (fig. 5(c)), and the induced rolling moment decreases by O.OlCx, (fig. 
5(e)) . The rapid turnover that occurs near y/b = 1 is a limitation of the lat- 
tice approach at this time, and is discussed later under "Effect of Wake Roll- 
Up" . 


This formation is roll stable, i.e., with flaperon fixed in a trimmed 
position, a spanwise movement of wing B would cause changes m induced rolling 
moment that would tend to return the wing to the trimmed conditon. The results 
for the spanwise scan with wing B at x/b - -3 are shown in figure 6. The in- 
duced drag ratio for wing A (fig. 6(a)) shows a rapid variation from the thrust 
value seen m figure 4(a); as y/b increases from 1.0 to 1.05 the ratio increas- 
es by 0.46. The kin k at y/b = 1.0 seen m figure 4(a) for wing B is present 
here also for wing A. 

The incidence and induced lift ratios for wing A (figs. 6(b) and (c)') also 
have rapid variation near y/b - 1.0; as wing B moves from y/b ™ 1.0 to 1.05, 
the incidence ratio for wing A increases by 0.07 and the corresponding induced 
lift decreases by .068 C L 

The rolling moment on wing B is not plotted for this case; the maximum 
flaperon deflection calculated at x/b = -3 is about 0.1°. 


Vertical Scans 


The results for a vertical scan at x/b - 3.0, y/b - 0.89, are given in 
figure 7. These show a strong sensitivity to vertical movement. The 0.21 in- 
duced drag ratio for wing B (fig. 4(a)) increases to 0.35 for a vertical move- 
ment of 0.05b, see figure 7(a). The incidence ratio, meanwhile, increases by 
0,04 (fig. 7(b)) and the induced lift ratio, AC^C^, decreases by 0.025 (fig. 
7(c)). 


The formations considered here keep the following wings (i.e., wing B) 
above the wake of the leader. This is stable as far as lift is concerned, 
ie.e, a vertical displacement upwards would reduce the induced lift, (fig. 
7(c)), causing the wing to come back down, and vice-versa. With wing B below 
the wake, a downward displacement would reduce the induced lift and so the 
wing would continue going down. 
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WING A 



FORMATION 

TYPE 



(a) Induced Drag Ratio. 


Figure 5. Effect of Spanwxse Position: x/b = 3; 

z/b = 0.0. 
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WING A 



(b) Incidence Ratio. 



(c) Lift Increment in Formation. 


Figure 5. Continued. 
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VERTICAL POSITION z/b 

(a) Induced Drag Ratio. 



(b) Incidence Ratio. 


(c) Lift increment in Formation. 


Figure 7. Effect of Vertical Position: x/b = 3; y/b = .89. 
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The effect of vertical position on flaperon deflection is shown in figure 
7(d) • A vertical movement of 0*05b from z ~ 0 causes a reduction in rolling 
moment from -.0134C L to -.00$5Cx,, (fig* 7(e)). In fact, the induced rolling 
moment changes sign (goes positive) above a z/b value of about 0.2 before re- 
turning towards zero. This is because the position of maximum vertical in- 
duced velocity on wing B moves outboard as wing B moves vertically; whereas 
the inner tip of wing B receives the maximum upwash on the wing when z/b - 0, 
the upwash there rapidly goes toward zero as the tip moves above wing A* f s trail- 
ing vortex sheet. 

The results for the vertical scan at x/b = -3 are shown in figure 8. In 
this case, the results are shown for negative z/b values, i.e., wing A is in 
the stable lift position above the wakes from wings B. Again, rapid changes 
with vertical displacement are indicated; a 0.05b vertical movement from z = 0 
causes an increase of 0.38 in the induced drag ratio on wing A (fig. 8(a)) 
while the incidence ratio increases by 0.05 (fig. 8(b)) and the induced lift 
increment decreases by 0.45C L , (fig. 8(c)) . The rolling moment for wing B is 
not presented for this case since it is constant and very small (£_ - 0.1°) . 


Effect of Wake Roll-up 


In the calculations discussed so far in this Section, the trailing vortex 
sheet shed by the leading wing (or wings in the case of negative x/b values) is 
represented by a parallel lattice of semi- infinite vortices. In real flow, the 
sheet would be almost completely rolled up at three spans downstream, the vor- 
tex cores being at the centroid of shed vorticity for each semispan. In the 
calculations where the lattices overlap (i.e., y/b < 1, see fig. 1), the trail- 
ing vortices that pass over the wing induce a downwash over the - tip region, 
whereas they would induce an upwash there if they were at the rolled-up 
vortex position. Neglecting roll-up, therefore, gives reductions m calculated 
induced lift and induced thrust; these effects are indicated in the calculated 
spanwise distributions of lift and induced drag in figure 9(a) and (b) , res- 
pectively. The calculations discussed so far are, therefore, pessimistic for 
the maximum induced drag saving, but larger rolling moments would be expected 
to occur with wake roll-up accounted for. 

The loss in induced upwash as the lattices overlap causes the kinks ob- 
served m the spanwise scan results at y/b =1.0 (figs. 5 and 6) . Also the 
turnover in the curves of incidence ratio, induced lift, flaperon angle and 
induced rolling moment (figs. 5(b), (c) , (d) , (e) and 6(b) and (c)) would occur 
at a smaller y/b value as the wing tip passed through the vortex core position. 

Two attempts were made to try to allow for the effects of roll-up in the 
formation with x/b = 3.0 and y/b = 0.89. The first attempt used the iterative 
vortex roll-up capability of the present calculation method. Several computer 
runs were carried out, but the calculations were not successful because of the 
close interaction between the opposing tip vortices and the wing tip. After 
two iterations, the • calculated vortex trajectories passed too close to the 
following wing control points for reliable force calculations to be made. 
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VERTICAL POSITION z/b 

(a) induced Drag Ratio. 


Figure 8. Effect of Vertical Position: v/h 






(b) Spanwise Distribution of Induced Drag. 

Figure 9. Spanwise Distributions of Lift and Induced Drag: 
x/b = 3; y/b - .89; z/b = 0. 


29 





The problem is essentially a close vortex/wing interaction, and this poses 
numerical difficulties for a vortex-lattice method. Recent developments in 
near- field modeling (refs- 7 and 8), would remove these difficulties, but the 
developments are not yet incorporated in the method. 

In a second attempt to allow for wake roll-up, the leading wing and its 
trailing vortex lattice were replaced by a single vortex having the total cir- 
culation of all the trailing vortices on the leading wing semispan. The vortex 
was placed 0.1b away from the following wing tip. This was done to try to 
avoid numerical difficulties that occur when calculating the induced thrust on 
a wing in the presence of a vortex using equation (2) and the vortex- lattice 
method. Even so, an unduly large induced thrust ratio, -0.62, was calculated. 

Further work is clearly required to calculate the effect of roll-up of 
these formation calculations. 3uch work should consider near-field modeling 
techniques such as developed in references 7 and 8. 
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OPTIMUM FORMATIONS 


The present calculations indicate that formation flying with all aircraft 
having the same induced drag saving is not practical because of the small dis- 
tances separating the wings. However, comparable savings for the formation as 
a whole are obtained by having streamwise separation of three or more spans, 
(fig. 4(a)). This would allow very open formation and should lead to safer 
operation. With large streamwise spacing, the leading aircraft is essentially 
m free air, so each aircraft in the formation would take a turn in the lead 
position to equalize fuel used. 

The reason for the improvement in benefit for the following wing as the 
spacing increases is that the upwash induced by the trailing vortices from 
the leading wing increases with downstream distance and quickly achieves double 
the initial value. (In effect, it is the transition from a semi-infinite vor- 
tex to an infinite one, and is demonstrated for a single vortex m figure 10.) 
Wake roll-up would enhance this effect. 

Two basic types of open formation appear feasible. The first type is 
similar to the "V" vormation used by birds, except larger streamwise spacmgs 
would be used. The formation could start in echelon (fig. 11(a)), then the 
leading aircraft could drop back after a time, and so on, until the echelon 
was reversed. Each aircraft would then have been in the lead position. All 
the aircraft would require trimming in roll in this formation. 

The second type of formation is based on the negative x/b positions, and 
would have two rows of aircraft. The rows could be swept (fig. 11(b)) to in- 
crease the separation between the aircraft m each row. Aircraft in the lead- 
ing row would have virtually no drag saving, but the second row would achieve 
double the saving. This type of formation has no basic roll- trim problem, and 
could be used for ferrying small aircraft. Large, long range aircraft could 
form the first rows, and each "slot" in between could support several smaller 
aircraft (e.g., R.P.V.'s) in line astern, (fig. 11(c)). Positions m each 
"slot" are roll-stable; they are also altitude stable if the following wings 
are above the level of the trailing vortices. 

For the echelon or "V" type of formation, the induced drag ratio for N 
aircraft would be 



C 


D 


0.2 


+ 


0.8 

N 


( 8 ) 
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(One leader with no saving and N-l followers with an induced drag ratio of 

0 * 2 .) 

For the double row type, the induced drag ratio for N (odd) aircraft would 
be 




0.26 


0.74 

N 


( 9 ) 


( (N+l) leaders with no saving and (N-1) followers with a thrust ratio of 
2 2 

0 . 48 .) 

The induced drag ratio for the double row formation could be reduced by 
adding an aircraft to each end of the second row. These aircraft would have 
the 0*2 induced drag-ratio, and would require the trimming in roll. 
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RANGE EXTENSION 


Any general conclusions relative to increase in cruising range which may 
be expected from reduced induced drag must be approached with caution because 
factors such as buffet boundaries and drag rise with Mach number are aircraft- 
dependent, Further, because of lack of simulation and flight experience to 
assess the feasibility of accurate stationkeeping in formation, it is not known 
what fraction of the potential reduction in induced drag is achievable under 
operational conditions. However, automatic stationkeeping using the aircraft 
automatic flight control system and laser or microwave position determination 
may be necessary, especially for flights of long duration. 

In view of this, it should be recognized that the estimated increases in 
cruising range using the relationships derived in Appendix A are potential 
increases, and must be scaled down by some factor depending on the character- 
istics or constraints of the aircraft considered. 


In estimating potential increases in range by flying in formation the fol- 
lowing assumptions were made: 


* 


* 


* 

* 


/ 

* 

* 


Cruise Mach number is the same in formation flight as for the 
case of a single aircraft not flying in formation. 

Induced drag can be approximated by the expression 
C D. " (S/ C L 2 ) C L 2 ' Wher * ( C d/ C L 2 ) iS a COMtant - 
Specific fuel consumption is constant. 

Altitude is allowed to change in order to achieve the desired 
values of 

C L /C D or 0 L‘ /C D- 


Mach number effects and buffet boundaries are neglected. 

Zero-lift drag, C , is constant and the same for each aircraft, 
D 0 


Aircraft in formation change position frequently enough that an 
average induced drag coefficient, C D ^, can be used for the forma- 
tion, and each aircraft will consume the same amount of fuel 
during the mission. 


With these assumptions, four formation flight 
and estimates were made for range of the formation 


strategies were considered, 
relative to a range of a 
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single aircraft not flying formation. These four strategies, and the corre- 
sponding single-aircraft strategies used for determining relative ranges, are 
as follows: 


A Formation at its max C_/C_; single aircraft at its max C /C . 

Xi D D 

h h 

B Formation at its max /C^; single aircraft at its max C^/C^. 


C Formation flying flight profile similar to that of a single 
aircraft flying at its max C^/C^. 


D Formation flying flight profile similar to that of a single 

aircraft flying at its max C^/C . 

L D 


Equations derived m Appendix A for each of the above cases were used to 
estimate the relative ranges for each of the above cases, and plotted in the 
figure below for an echelon formation; double- row formations yield very nearly 
the same results - 

In Figure 12, it can be seen that greatest range increases occur 
when the formation flies flight profiles optimized for the formation (cases A 
and B) rather than when the formation flies flight profiles similar to those 
optimized for an aircraft not flying m formation (cases C and D) . 

The figure shows significant befit can be derived even for formations of 
only two aircraft (as much as 30% increase for cases A and B) . A point of 
diminishing return is reached beyond formations of ten aircraft* A reasonable 
trade-off between a practical formation size and range benefit seems to lie at 
about three to five aircraft, with corresponding maximum potential range in- 
creases of about 46% to 67% (cases A and B) . Corresponding increases for case 
D for three to five aircraft formations are about 15% to 19%. Even these in- 
creases could be significant for strategic augmentation missions. Very large 
formations, say 15 aircraft, show maximum potential increases of about 100%. 

Again, a word of caution: it is not known what fraction of this potential 
range increase is achievable in practice. 


36 



2 


3 4 5 


IO 


20 30 40 50 


NUMBER OF AIRCRAFT IN FORMATION 


Figure 12. RELATIVE RANGE VS. NUMBER OF AIRCRAFT IN FORMATION 
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CONCLUSIONS AND RECOMMENDATIONS 


Conclusions 


1) Formations giving equal induced drag saving on all aircraft require 
dangerously small spacings, but comparable savings for the formation as a 
whole are obtained with large streamwise spacing. Leading aircraft in 
the open formations have virtually no induced drag saving. 

2) Two types of open formations are possible. The first type is a "V" or 
echelon type having one leader. Each of the following aircraft would 
require trimming in roll. The second type of formation has two rows of 
aircraft. The leading row has virtually no induced drag benefit, but 
the second row has double the benefit of the echelon type. The second 
type of formation has no basic out-of-balance rolling moment. 

3) Significant induced drag savings are acheived in both formations. Air- 
craft in the echelon formation have an induced drag of the order of 20% 
of the free- air value at the same C L . The second row of the double row 
formation has an induced thrust of 48% of the free- air induced drag at 
the same C^. 

4) Both formations are insensitive to longitudinal movements, but benefits 
(and trimming problems) decrease rapidly with lateral or vertical move- 
ments from the optimum position. 

5) A reasonable tradeoff between a practical formation size and range benefit 
seems to lie at about three to five aircraft with corresponding maximum 
potential range increases of about 46% y to 67% with a flight profile opti- 
mized for the formation (rather than one optimized for aircraft flying 
alone) . At this time it is not known what fraction of this potential 
range increase is achievable in practice. 


Recommendations 


Maximum benefits are obtained when the vortex core from the leader passes 
close to the tip of the follower. This is a close vortex/wing interaction 
problem which is beyond the scope of present calculation methods. It is recom- 
mended that the numerical problems of calculating close vortex/wing interaction • 
be investigated, based on recent developments in near-field techniques (refs. 

7 and 8) . Such an investigation should lead to a more accurate assessment of 
the maximum benefits, and of the sensing and control aspects near the optimum 
locations . 


38 



Additional calculations should be done for different lattice densities. 

The computer program size limited the number of vortex quadrilaterals for the 
present calculations. Extensions of the program to allow more quadrilaterals 
should be considered to investigate the convergence of the results. The method 
should then be used to assess the range benefits and stability and control 
aspects of a specific aircraft in formation. 
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APPENDIX A 


DERIVATION OF CRUISING RANGE EQUATION FOR A JET AIRCRAFT 


NOTATION : 

S = wing area 

R = aircraft cruising range, miles 
V = true airspeed, MPH 

C <* specific fuel consumption, lbs, per hour per pound of thrust 
L = lift 
D = drag 

W = aircraft weight 

W = aircraft weight at beginning of cruise 

D 

W = aircraft weight at end of cruise 

ft 

a = air density at altitude divided by sea level density 
p - air density 


( )_ = Denotes formation of aircraft 

Jb 

( )pgL = Denotes following aircraft 
( ) = Denotes lead aircraft 


( )g = Denotes a single aircraft (flying alone) 



The purpose of this appendix is to develop some measure of potential 
range increase which may result from reduction in induced drag due to flying 
in formation. 

It is recognized that the benefits of flying in formation are very 
much aircraft dependent, and that constraints such as buffet boundaries 
and drag rise as a function of Mach number together with control problems 
will probably prevent achieving full potential benefit. However, some 
generalized computations would still be useful to provide some insight for 
establishing range performance goals. 


A. 1 FORMATION CRUISE AT ITS MAXIMUM LIFT-DRAG RATIO 

Consider first the case of a single aircraft not flying in formation. 
Assume that it normally cruises, for maximum range, at a Mach number below 
the drag rise, perhaps Mach 0.6 to 0-7, and at or near its maximum lift-drag 
ratio. With these assumptions the cruise altitude will be a function of 
wing loading (or gross weight, since wing area is fixed). 


The classical range equation for jet aircraft is given by the follow- 

ing equation: R = “f E tt 7 ^ % (A-l), 

Jw B C C D W 


Assume V, C, and C ,/C are constant, then 
C L D 

K ■ c ^ 108 <vv 


(A-2) 


Now, define relative range as the maximum range of a formation of air- 
craft divided by the maximum range of a single aircraft. Assuming changes 


A-2 




In cruise Mach number, specific fuel consumption and zero-lift drag are 
negligible, the range of a formation of aircraft at its optimum cruise 
altitude relative to that of a single aircraft at its optimum cruise 
altitude becomes: 





(A-3) 


The task now is to determine C T /C for the single aircraft and for the 

t» D 

formation. Assuming that the equation for drag coefficient takes the form 
C D = C Do + C Di 

= C D 0 + < C Di/ G L 2 >C L 2 • (A-4) 


it can be shown that C^/C^ is maximum when Cj^ = C])q = (Cp^/C^ )C^ 


The drag coefficient for this condition is then written 

S - 2 % - 2 «V c l 2 > c i . 2 


And the corresponding lift coefficient is 

c ^ 

L 

Subsituting 

‘Wnux ' <A - 5) 


Using the subscripts S and F to distinguish the lift-drag ratio for a 

single aircraft not flying in formation, and for the average of a number 

of aircraft in formation, respectively, the relative range equation becomes, 

after substitution for (C T /C ) , 

v L D max 
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(A- 6 ) 


*F < C V C I.V 

R s (V^V 


Assuming that the lead aircraft in the formation derive no change in 
induced drag from the follower aircraft, i.e., their aerodynamic characteris- 
tics are basically the same as if they were flying as single aircraft, and 
assuming the formation may contain' lead aircraft, follower aircraft, 

and N total aircraft (lead plus follower aircraft), the average induced drag 
coefficient for the formation is given by. 


N L C Di + n FOL C Dj 


(A-7) 


Since the entire formation will fly at the same lift coefficient, and 
since the drag coefficient of the lead aircraft will be the same as for a 


single aircraft not flying in formation, 

(C n ,/C 2 ) - N L (C Di /C L + N F 0 L (Cd 1 /C L ^FOL 

1 L r „ 


(A-8) 


Substituting into Equation A-6, 


R s L 


N L + N FOL (C Di /C L 2) FOL /(C Di /C L 2) S 


"1 

1 


(A-9) 


C ^s/C 

A. 2 FORMATION CRUISE AT ITS MAXIMUM L 2 ' D 


Recall the range equation (Equation A-l) 


J V dw 

c C D w 


In equilibrium flight 

fo T.T 1% 


*2 W T 

P 0 <JSC L 


(A-10) 



Assume specific fuel consumption, C, is constant. Substituting for V in 


the range equation 




R = 


1 ( , 2W M dw 

“ C J SO ' 1 c 


C D W 


(A-ll) 


If drag can be approximated by the expression = 

l, 

the ratio (C , vC ) is maximum when 
L 0 _J< 

% 


C L" 


3(C Di ) /c l 2 


% + \ ,c lW- 


(A-12) 


The corresponding drag coefficient can then be written 

V 

C = (4/3)C (A-13) 

u u Q 

and, combining equations A-12 and A-13, 

“vVmax - (1/4) « /C D < C d. /C L 2 >' , ‘ (A - 14) 

0 l 


Assume that as fuel burns the aircraft is allowed to climb at constant 

h 

angle of attack, i.e., C and (C ^)/C_ remain constant, and that cruise Mach 

L L U 

number and velocity remain constant. This being the case W/ais constant and 
1 S W , 

R ■ C < C lX> l0S <VV . <A ~ 15) 


Range of a formation of aircraft relative to a single aircraft then 
becomes 


V R S ■ < V V^< c lXv (c lXV 


(A-16). 


is determined from the relationship 


W = ho Vo C 
W B ^ ° B L 


CA-17) 
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If V is essentially the same in going from single aircraft flight to 


formation flight 


(o bVf “ ^bVs 


and 


/° R - C t /C L 

B s b f l f l s 


Recall that for maximum CL /C 

'■ C Dn ■ 


i«V c L 2 j 


(A-18) 


Then 


g B S = (C Dj /C L 2) S^ 
c bf " (c Di /c L 2 ) F ^ 


(A-19) 


Substituting expressions for C^/C^ and into the relative range 

(Equation (A-16), 


■ *F 

R S (C^V 


(A-20) 


^2 # 

Interestingly, the relative range obtained by maximizing /C^ is the 
same as that obtained by flying maximum C^/C^, and the relative range for the 


formation can again be written 


h p L + N : 

R o - 


FOL (CD i /C L 2) FOL /(C Di /C L 2 


N 




-h 


(A-21) 


A. 3 FORMATION CRUISE PROFILE SIMILAR TO THAT OF A SINGLE AIRCRAFT FLYING AT 
ITS MAXIMUM LIFT DRAG RATIO 


Assume that a single aircraft not flying in formation flies at a maxi- 
mum, constant lift-drag ratio, and that the formation flies a flight profile 
similar to that of a single aircraft not flying in formation. This means 


A-6 



that the range of the formation relative to that of an aircraft flying single 
at its maximum lift-drag ratio can be expressed as follows: 


b 

R S 

b 


- w 

“ ( V V ( vv 

1 + % 'V 

_____ s o 

1 + 


CD i F /C D f 


(A-22) 


For maximum lift-drag ratio (single aircraft). 
Hence, 


V*S 


1 + < C V°V 




(A-23) 


2 2 

Assuming that Cj^ can be approximated by )C^ and that, 

because the flight profile of the formation is assumed to be similar to that- 
of a single aircraft, the lift coefficients are the same, then 


R s 1 + (0 Dl /c L 2 ) F /(c Dl /c L 2 ) s 


_(A-24) 


A. 4 FORMATION CRUISE PROFILE SIMILAR TO THAT OF A SINGLE AIRCRAFT FLYING AT 
ITS MAXIMUM C T h / Cp ~ 


By analogue to the case of a formation cruise profile similar to that 

of a single aircraft programmed to fly at a constant, maximum C^/C^, relative 

% 

range for flight at maximum C /CL can also be written 

L D 


Rp 1 + ° D i„ /C D 

rT 


0 


Cn /C 
F 0 


1 + ”D 


(A-25) 
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k, 

From Equation A-13, for maximum C /C , 

L> D 

C D = 4/3 C D 0 


Writing 

C D ±S = 1/3 C D 0 

(a) 


c Dq = 3 ^ D i s 

(b) 


C D ±F ■ (c Dl /c L 2 ) F c L 2 

(c) 


C„i s » (C Di /C L 2 ) s Cj 

(d) 

and substituting into Equation A-24 



4/3 


R s“ 1 

+ (l/3)(C Di /C L 2) F /(C Di /C L 2) 

S 


(A-26) 


(A-27) 
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RECAP 


The equations derived in this appendix are summarized below for con- 
venience and comparison. Assumptions made in deriving them are: 

Specific fuel consumption is constant. 

Velocity and Macli number are constant and the same for formation and 
single aircraft flight. 

Altitude is allowed to change to achieve the desired value of 

C lX ° r C I>D- 

Change in sound velocity with altitude is neglected. 

Mach number effects and buffet boundaries are neglected. 

2 2 

Induced drag can be approximated by the expression Crw =* (Crj./C )C T , 

u l i» h 

2 

where (Cjj^/C^ ) is a constant. 

Profile drag, C-q^, is constant and the same for each aircraft. 

Induced drag of the lead aircraft of the formation is not affected by 
other aircraft in the formation; follower aircraft do experience changes in 
induced drag from aircraft immediately preceding. 



N 

AVERAGE INDUCED DRAG COEFFICIENT FOR A FORMATION OF AIRCRAFT WITH L LEAD- 
ERS AND N F qL FOLLOWERS 


g D-! = N L Cp iL CD iFQL 


N 


, p 2. _ N L <C »1 /0 l\ + »FOI. < C »i /C L 2 >FOL 

(C Di /C L >F N 


(A-28) 


Assuming 

< c d^ c l 2 >f 

< c V C L 2 >S 


' (C I>i /C L 2) C L Z > and C D± l * CD i s - 
H L + %L (C °i /C L 2 >F0L /( % /C l. 2:i S 


N 


then 


(A-29) 


RANGE OF A FORMATION CRUISING AT ITS MAXIMUM FORMATION LIFT-DRAG RATIO 
RELATIVE TO THAT OF A SINGLE AIRCRAFT FLYING AT ITS MAXIMUM LIFT-DRAG 
RATIO ~ ' 


h <-W c l\ h 

R s <CD 1 /C L 2 jj J !:i 


(A-30) 


RAMGE OF A FORMATION CRUISING AT ITS MAXIMUM Ci.^/Cn RELATIVE TO THAT OF A 
SINGLE AIRCRAFT FLYING AT ITS MAXIMUM C T ,%/Cn 


h <°D±/<%h k 

R s (c Dl /c L \ h 


(A- 31) 


RANGE OF A FORMATION FLYING A CRUISE PROFILE SIMILAR TO THAT OF A SINGLE 
AIRCRAFT FLYING AT ITS MAX CT./Cn 


R S ~ 1 + ( C Di/ C L 2) F /(C Di /C L 2) S 


(A-32) 


RANGE OF A FORMATION FLYING A FLIGHT PROFILE SIMILAR TO THAT OF A SINGLE 
AIRCRAFT FLYING AT ITS MAX C T h/Cn 


% _ 4/3 

R g 1 + (l/3)(C Di /C L 2) F /(C DjL /C L 2) s 
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NUMERICAL EXAMPLES 


Representative values obtained from Section 5 for use in estimating 
relative range potential are as follows: 



FORMATION ' 

ECHELON 

DOUBLE ROW 

° D W CD is 

0.2 

-0.48 


Substitution of these values into Equations A-27 through A-32 yields 
the results summarized in Table A-l. Some comments on these results are as 
follows : 

o Echelon formations provide slightly greater benefits than 
double row formations. 


o Relative ranges are significantly greater for formations 
cruising at flight conditions optimized for maximizing 
range of the formation (cases A and B in Table A-l) rather 
then flying flight profiles optimized for maximizing range 
of a single aircraft. 

o Potential cruising range for a formation of five aircraft in 
estimated to be 67 per cent greater than that of a single 
aircraft. Constraints, such as buffet boundaries, and 
possible difficulties in stationkeeping control would 
probably reduce this potential benefit substantially. 


A-ll 






N 




r f /r s 




ECHELON FORMATION 

DOUBLE ROW FORMATION 

E19 

A & B 

C 

D 

Ell 

A & B 

: 

C 

D 

CD4 

c Di 

1 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

2 

.6 

1.29 

1.25 

1.11 

.63 

1.26 

1.23 

1.10 

3 

.467 

1.46 

1.36 

1.15 

.507 

1.4 

1.33 

1.14 

5 

.36 

1.67 

1.47 

1.19 

.408 

1.57 

1.42 

1.17 

7 

.314 

1.78 

1.52 

r. 2 i 

.366 

1.65 

1.46 

1.19 

9 

.298 

1.86 

1.55 

1.22 

.342 

1.71 

1.49 

1.20 

15 

.253 

1.99 

1.60 

1.23 

.309 

1.80 

1.53 

1.21 

39 

.221 

2.13 

1.64 

1.25 

.279 

1.89 

1.56 

1.22 

99 

.208 

2.19 

1.66 

1.25 

.267 

1.94 

1.58 

1.22 

INFINITE 

.20 

2.24 

1.67 

1.25 

.26 

1.96 

1.59 

1.23 


TABLE A-l . RELATIVE RANGE AS A FUNCTION OF NUMBER 
OF AIRCRAFT IN FORMATION 

A = Formation cruising at max formation C|_/Cp 
Single aircraft cruising at its max fy/Cp. 

B = Formation cruising at max formation C^/Cq. 

Single aircraft cruising at its max C^/Cp. 

C = Single aircraft cruising at its max C|_/Cp. 

Formation cruising at single aircraft's max Ci/Cp 

D = Single aircraft cruising at its max Cl 2 /Cd. 

Formation cruising at single aircraft's C|_ 2 /Cp. 
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